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Abstract - Absorbing phase transition in restricted exclusion processes are characterized by 
simple integer exponents. We show that this critical behaviour flows to the directed percolation 
(DP) universality class when particle conservation is broken suitably. The same transition, when 
studied using the average density as the controlling parameter, yields critical exponents quite 
different from DP; we argue that these exponents are actually related to DP by a scaling factor 
1//3dp- These conclusions also apply to conserved lattice gas in one dimension. 



Directed percolation(DP) [I] is the most dominant uni- 
versality class of absorbing phase transition (APT). The 
famous 'DP conjecture' [2] claims that that any APT with 
a fluctuating scalar order parameter should generically 
belong to DP class. Several models systems, like reac- 
tion diffusion processes [3] , depinning transitions [4] , dam- 
age spreading [5] , synchronization transition [6] , sand-pile 
models [7] and certain probabilistic cellular automata [8] 
are known to undergo APT belonging to this universality 
class. Recently the DP critical behaviour has been verified 
experimentally in context of liquid crystals [9] . 

Surprisingly, in spite of its robustness directed perco- 
lation does not have an exact solution, even in (1 + 1)- 
dimension, although the critical exponents are calculated 
to high accuracy JTO] . Only a few models of APT are 
solved exactly; examples include compact directed per- 
colation (CDP) branching-annihilating random walk 
(BARW) [12] , conserved lattice gas (CLG) model in one 
dimension [13] and restricted asymmetric simple exclu- 
sion process (RASEP) [14]. The non-DP behaviours in 
CDP and BARW owe to the special symmetries of the 
order parameter, namely particle-hole symmetry and par- 
ity respectively. The scalar order parameters of CLG and 
RASEP do not have any special symmetry but it is cou- 
pled to an additional conserved field present in the system. 
Conventionally, presence of this additional field is blamed 
for these models to be in a non-DP class. This reasoning 
is consistent with the non-DP behaviour observed in fixed 




Fig. 1: Propagation of activity in asymmetric (left) and sym- 
metric (right) restricted exclusion process starting from the 
maximally active configuration {110110...}, for w = 0.560. 
The blue (dark) pixels represent active sites. 

energy sand pile models [15] and a conjecture by Rossi et. 
al. [16] that critical behaviour of any APT in presence of 
a conserved field should differ from DP. 

However, more recent investigations reveal [17] that the 
actual critical behaviour of one dimensional fixed energy 
sand-pile models belong to DP. The apparent non-DP be- 
haviour was only a transient phenomenon; coupling of 
the order parameter to the additional conserved height 
field does not necessarily induce a new critical behaviour. 
APTs in energy exchange models [18] too show DP be- 
haviour even though the total energy of the system is con- 
served. In view of these examples it is not convincing to 
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accept 'conservation' as the cause of an absorbing transi- 
tion not being in DP, unless one explicitly shows that the 
specific model flows to DP if the existing conservation is 
broken. 

In this article we study the restricted asymmetric sim- 
ple exclusion process (RASEP) as an example of a simple 
and exactly solvable non-DP critical behaviour in pres- 
ence of a conserved density field. The dynamics of RASEP 
is extended to incorporate deposition and evaporation of 
particles in the bulk so that the conservation of density 
is violated; the system now undergoes an APT at a crit- 
ical deposition rate w c . The resulting critical behaviour 
is found to be in the DP universality class. However, the 
critical exponents turn out to be much different from DP 
when the average density p is treated as the control pa- 
rameter, like the conserved RASEP. We show that these 
exponents are actually related to DP by a scaling factor 
1/ ' Pdp as the singular behaviour of p at the critical point 
is also characterized by DP exponents. 

Model. — The restricted asymmetric exclusion pro- 
cess is defined on a one dimensional periodic lattice la- 
beled by sites i — 1,2... L. The sites are either vacant or 
occupied by at most one particle; corresponding site vari- 
ables are Sj = 0, 1. A particle from a randomly chosen site 
hops forward to its neighbouring vacant site only when it 
is preceded by at least p particles. Thus, the model with a 
specific p evolves according to the corresponding particle 
conserving dynamics, 

p = l : 110 -S- 101 
p = 2 : 1110^1101 

p = 3 : 11110^11101 ... etc. (1) 

Any configuration which does not have more than p con- 
secutive particles anywhere in the system is an absorbing 
configuration. If the system reaches such a configuration 
it gets stuck there forever. All the other configurations, 
which have at least one block of more than p consecu- 
tive particles, are active. Since all the configurations with 
particle density p > p/(l + p) are active, the APT is ex- 
pected to occur at some lower density. It has been shown 
analytically [14] that the system undergoes an APT at 
p c = p/(l + p) exactly with critical exponents 

/? = 1,^ = 1,77 = 1,/?' = 0. (2) 

Moreover, numerical simulations have suggested |19| that 
a = 1/2, z = 1, and v\\ = 2. 

To claim that the conserved density is the sole cause of 
the non-DP behaviour of RASEP (Eq. ©) we intend to 
break the conservation of density and investigate whether 
this truly drives the transition to the DP fixed point. How- 
ever, any arbitrary dynamics that deposits (evaporates) 
particle at (from) a vacant (occupied) site may altogether 
wipe out the possibility of a phase transition of the sys- 
tem. For example, the simplest way one can break the 
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Fig. 2: Decay of p a (t) from a maximally active configuration 
{110110110...} for different values of w. The same data is 
plotted as p a {i)t a versus t in the inset. In the long time, this 
curve becomes horizontal for w c = 0.56757 . Here, L = 2 16 . 

density conservation is by augmenting the dynamics with 

l-w 

Otil, 

i<- 

This bulk dynamics fixes the average density of the sys- 
tem p = w. However, there are no absorbing configura- 
tions left anymore as the dynamics allows any arbitrary 
configuration to become active eventually . Thus, the sys- 
tem is active for any w and does not show any transition. 
To preserve the APT, the non-conserving dynamics must 
be designed suitably, so that the absorbing states are still 
present in the phase space. One possible way to ensure 
this is to device a dynamics which does not affect the ex- 
isting absorbing configurations. In other words, one ought 
to deposit or evaporate particles only on active configura- 
tions - all the configurations which were absorbing in the 
conserved model would then still be absorbing in the non- 
conserved version. 

A possible non-conserving local dynamics for RASEP, 
which fulfills this condition for any p, is 

1 — w 

p = 1 : 110 111, 

w 

l-w 

p = 2 : 1110 1111, 

W 

1—w 

p = 3 : 11110 11111 ... etc, (3) 

w 

where < w < 1 fixes the density by controlling both the 
deposition and evaporation of particles. Clearly, large w 
would correspond to high density. This dynamics, when 
added to the conserving evolution rule ([l} , allows addition 
or removal of a particle in the neighbourhood of an active 
site only. Absorbing configurations, which do not have 
(/i + 1) or more consecutive particles, are not affected by 
the non-conserving dynamics (|3|). 

Although the number of particles in the system fluc- 
tuates, the average density p = (s^ is a non-decreasing 
function of the deposition rate w. For small values of w, 
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Fig. 3: The critical exponent /3: (a) Plot of p a (t) vs. t for 
different values of w - w c = 0.0001, 0.0002, .., 0.0512 for a sys- 
tem of size L = 2 15 . (b) Corresponding saturation values are 
plotted against w — w c in log-scale. The best fit straight line 
corresponds to the slope /3 = 0.278(1). 

the system is likely to get arrested in some absorbing con- 
figuration with number of particles N < -^pj. Again, for 
w ~ 1, the density fluctuates about a large average value 
p(w) (note that p(l) = 1) and the system remains active. 
Thus, we expect a transition from the active phase to the 
absorbing state when w is decreased below a critical value 
w c , such that p{w c ) < T^pi- 

Let us discuss the \x = 1 case in details. This exer- 
cise can be repeated straightforwardly for any p > 1. The 
complete dynamics for RASEP with p = 1 is 



110 -> 101 



conserving 



110 



1 — W 



111 non — conserving. 



(4) 
(5) 



Any configuration devoid of 110 is absorbing. In other 
words all the configurations which have two or more con- 
secutive Is are active. Thus, 



Pa = (Si-XSi(l - Sj+l)) 



(6) 



plays the role of order parameter. In the following we use 
Monte-Carlo simulations to study the critical behaviour of 
this APT. 

Critical point and decay exponent a: Starting from 
a homogeneously active configuration, the activity p a [t) ~ 
t~ a at the critical point w c . The propagation of activity 
is shown in the left panel of Fig. [TJ In the active phase 
w > w c , the activity eventually saturates to some finite 
value and in the absorbing phase p a (t) vanishes in the long 
time limit. One can then estimate the critical point w c 
and the exponent a by plotting p a (t) versus t for different 
values of w and looking for a power law decay. These 
estimates can be verified by plotting p a (t)t a against t for 
different values of w; the curve which remains horizontal 
in the long time limit corresponds to the critical point w c . 
This procedure, illustrated in Fig. [5J leads to the estimates 

to c = 0.56757(1) and a = 0.15984(7). 

The decay exponent a is in very good agreement with the 
corresponding DP value (see Tabic [TJ). The density at the 




Fig. 4: The data collapse according to Eq. |[8j for \w — w c \ = 
0.0016,0.0032,0.0064 using a = 0.1598 and v n = 1.74. The 
inset shows the unsealed data. The system size is L = 2 15 . 



critical point is 



p(w c ) = 0.4915(5), 



which is slightly lower than the critical density p c = 1/2 
for conserved RASEP (i.e. when dynamics (|5|) is absent). 
The detailed variation of p with w will be discussed later 
in this article. 

The order parameter exponent j3: In the super 
critical regime w > w Cl the activity saturates to some 
finite value p a - This saturation value shows a power law 
behaviour near the critical point, 



p a ~ (to - WcY 



(7) 



where /3 is the order parameter exponent. Figure EJa) 
shows plot of p a (t) versus t for various values of to — w c 
in the active phase. Corresponding saturation values are 
plotted against w — w c in Fig. EJb). The slope of the 
resulting straight line gives an estimate 



P = 0.278(1), 

consistent with Pop- 

Off-critical simulations and v u 

activity p a {t) obeys a scaling form 

Pa (t,w) =r a T{t/Q. 



For to ^ w c , the 



(8) 



where is the characteristic time scale of the system. 
Near the critical point, both in the active (w > to c ) and 
absorbing (to < w c ) phase this scale is expected to diverge 
as 



6, 



(9) 



The exponent v n can, thus, be estimated by plotting 
Pa{t)t a versus t\w — w c \ u " for different values of w and 
searching for a data collapse. This is shown in Fig. [4j 
the upper and lower branches correspond to w > w c and 
w < w c respectively This gives us an estimate 

v n = 1.74(1), 
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Fig. 5: Finite size scaling :p a (t)L^ l// "- L as a function of the 
scaled variable t/L z for systems of sizes L = 2 s , 2 9 , . . . 2 14 are 
collapsed using p/i/± = 0.254(4) and z = 1.58(1). The un- 
sealed data are shown in the inset. 

which is also in good agreement with the corresponding 
DP value 1.733847(6). 

Finite size scaling : For a finite sized system the 
decay of activity p a (t) at the critical point w c is expected 
to follow a scaling form 



Pa (t) = L-^g(t/L*), 



(10) 



where z is the dynamic exponent. Thus, at the critical 
point the plots of p a {t)L^/ u± against t/L z , for different 
values of L, are expected to collapse onto a single function. 
The best data collapse, as shown in Fig. [SJ is obtained by 
using 

P/vx. = 0.254(4); z = 1.58(1). 

The inset there shows the unsealed data for L = 2 8 -2 14 . 
Both these estimates are consistent with the DP expo- 
nents. It should be noted that this value of is 
strikingly different from that of the ordinary conserved 
RASEP, where j3/i/± — 1 was known exactly [T3] ■ 

The estimated exponents of non-conserved RASEP are 
summarized in Table [T] along with the corresponding expo- 
nents of DP for comparison. All the exponents are in excel- 
lent agreement with those of DP, indicating that RASEP 
indeed flows to the DP universality class when the density 
conservation is broken. 

In ordinary RASEP, particle density was varied to ob- 
serve an absorbing transition at p c = 1/2. It is therefore 
quite appropriate to ask whether the average density p, 
which now depends on w, can also be treated as control 
parameter in the non-conserved version. To answer this 
question, let us first look at the behaviour of p near the 
critical point. 

The density p(w). — The average particle density p 
is a well defined function of the deposition rate w in the 
supercritical phase w > w c . Near the critical point, how- 
ever, like all other observables the density p(w) is expected 
to have a non-analytic form, 

p(w) = p(w c ) + (w — w c ) K + higher order terms, (11) 
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:(a) Plot of p(w c ; L) versus l/L 



for L = 2 s — 2 15 . The inset shows plot of same date versus L -0 ' 6 ; 
linear fit yields p(w c ) = 0.4910. (b) p(t) is plotted against t in 
log-scale for various trial values of p(w c ); the best estimate 
p(wc) — 0.4915(5) results in a straight line. 

where p(w c ) is the value to which the density saturates 
at the critical point and the exponent k is related to the 
critical exponents, possibly in a non-trivial way. 

To estimate k we need to determine p(w c ) accurately. 
This can be done by extrapolating p(w c ; L), the saturation 
value of density for a finite system of size L, to the ther- 
modynamic limit. As illustrated in Fig. [6j p(w c ; L) ap- 
proaches a finite value in this limit l/L — Y 0. Since p(w c ; L) 
is not linear in l/L, it is useful to plot p(w c ; L) against L~ b 
with a suitable exponent b that gives a straight line. The 
inset of Fig. |6|b) shows that b — 0.6 does this job; a linear 
fit now gives an estimate p(w c ) = 0.4910. 

A more accurate estimate of p(w c ) can be obtained from 
the decay of p(w c ,t). Since, the characteristic timescale 
diverges at the critical point, in the long time limit we 
expect the reduced density p(w c ,t) = (p(w c ,t) — p(w c )) to 
decay algebraically 

pK,t)~r". 

Then p{w c ) can be determined by plotting p(w c ,t) ver- 
sus t in double logarithmic scale and adjusting the value 
of p(w c ) in a way that the resulting curve is a straight 
line. This is shown in Fig. Eta), and we have an estimate 
a = 0.1596(1) and p(w c ) = 0.4915(5). We proceed further 
with this estimate of p(w c ). 

Determination of the exponent k is described in Fig. 
[71 where saturation of p(w,t) for different values of w 
is shown in the left panel. Corresponding p~(w) obtained 
from this graph is plotted against w — w c in logarithmic 
scale in Fig. [JJb). The slope of this line over the first two 
decades corresponds to 

k = 0.278(3). 

Surprisingly, the exponents a and k, within the error bars, 
are identical to the DP exponents chop and /?dp respec- 
tively. This leads us to surmise that p(w) also undergoes 
an absorbing phase transition at w = w c with a decay ex- 
ponent a = a and order parameter exponent f3 = n. We 
investigate the other critical exponents to confirm that 
this transition, where the reduced density p(w) is the or- 
der parameter, also shows DP behavior. 
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Fig. 7: Determination of k : (a) The density p(w, t) versus t sat- 
urates to p(w) for different w - w c = 0.0001, 0.0002, .., 0.0512. 
(b) Corresponding reduced density p(w) — p(w) — p{w c ) van- 
ishes as (w - w c ) K with k = 0.278(3). Here L is 2 15 . 

First, the off critical scaling. In the inset of Fig. [5(a) 
we have plotted p{w,t) versus t for different values of w. 
Unlike w > w c regime where the saturation density p(w) 
varies with w, in the sub-critical phase p(w < w c ) = p(w c ), 
ensuring that the order parameter p(w) vanishes in this 
phase. The main figure Fig. HJa) shows the plot of 
p{w, t)t a versus t\w — w c \ v " for different values of w, which 
could be collapsed using 

i>„ = 1.73(1). 

This value of u n , again, is compatible with DP. 

At the critical point, p{w,t)L^/ u± - versus t/L z plot (see 
FigH^b)) for different values of L collapse following a scal- 
ing form similar to Eq. (fTOj) for 

P/ujl = 0.256(4); z = 1.58(1). 

All these exponents are listed in Table [TJ The fact that 
these exponents are consistent with those of p a as well 
as DP, leads us to conclude that in absence of particle 
conservation the reduced density p(w) can be regarded as 
an alternative order parameter. 

Density p as a control parameter. To have a 
direct correspondence with the APT in conserved RASEP, 
where density was tuned, we study the transition of the 
non-conserved model considering the average density p as 
the tuning parameter. Evidently, the critical point w = 
w c translates to p c = p(w c ) in the density axis. In the 
following we argue that the critical exponents of this p a 
versus p transition arc different from DP, but they can be 
expressed in terms of DP exponents in a trivial way. 

Along the density axis, the distance from the critical 
point is A = p—p c ; since p — p c ^ (w — w c )^, we have (w — 
w c ) ~ A 1 /' 3 . Therefore, critical exponents corresponding 
to those observables which vary algebraically with w — w c 
near the critical point (e.g. p a ~ A* 9 , ~ A~"" and 
~ A _!yj -) are now rescaled by a factor 1/(3dp, 



1; 



-IP- 



(12) 



Obviously, the decay exponent a 1 — P' /v' n and the dy- 
namic exponent z' = v'Jv'^_ do not change. 



Fig. 8: (a) p(w, t)t a for different values of w collapse when 
plotted against t\w — w c \ v " yielding an estimate v n = 1.73(1). 
(b) p(w c , t)L^I Vl - versus t/L z for different!/ are collapsed to the 
finite size scaling function for z — 1.58(1) and /3/i/± = 0.256(4). 
The insets of (a) and (b) show the unsealed pit) versus t data. 



The rescaling of exponents Eq. (|12|) can be confirmed 
by re-plotting the data of Fig. |3jb) and Fig. |4] using 
p{w) — p(w c ) instead of w — w c , which is shown in Fig. 
[9] The left panel of the figure shows p a as a function of 
p — p c resulting in the order parameter exponent /3' = 1. 
Data collapse according to scaling form Eq. (O is shown 
in the right panel for the supercritical regime p > p(w c ) 
with an unusually large v' n = {v n / P)dp = 6.269. Note that 
in this scenario, sub-critical branch does not exist as the 
density p does not vary for w < w c . 

It is worth mentioning that although the critical ex- 
ponents of this transition, when viewed by tuning p, are 
widely different from those of DP, the basic underlying 
critical behaviour is DP. It is only a coincidence that, like 
the conserved RASEP the order parameter p a here varies 
linearly with p — p c near the critical point. 

Restricted Symmetric Exclusion Process. Re- 
stricted Symmetric Exclusion Process (RSEP) is a varia- 
tion of RASEP where a particle is allowed to move in for- 
ward (backward) direction when it is preceded (followed) 
by at least p particles. In particular, the dynamics for 
p = 1, which is identical to that of 1-dimensional con- 
served lattice gas (CLG), is 



110 ->■ 101 <- 011. 



(13) 



It was shown [13l[T4] that this system also undergoes an 
APT with integer exponents ([2]) and belongs to the same 
universality class as RASEP. This non-DP behaviour seen 
in one dimensional CLG is very different from other re- 
lated models with an additional conserved field, namely 
conserved threshold transfer process (CTTP) and con- 
served Manna model though all three models in two and 
higher dimensions show identical critical behaviour. In 
ID, where CLG is exactly solvable, CTTP and Manna 
model shows a nontrivial absorbing phase transition be- 
longing to DP. It is natural to ask, what makes CLG dif- 
ferent from others ? 

One possible cause, like in RASEP, is the presence of the 
conserved density field. To confirm this we add particle 
deposition and evaporation 



110 



111 



011, 



(14) 
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Fig. 9: Exponents of p a vs. p transition: (a) An excellent 
match of p a vs. p — p c data (symbols) with a line of unit slope 
confirms /3' = 1. (b) Scaling collapse of p a (t) following Eq. {SJ, 
using v' K = (v„/P)dp- Inset of (b) shows unsealed p a (t). 



to the dynamics (fl3| of CLG. Clearly, this dynamics pre- 
serves the absorbing configurations of the system, assuring 
that there is still an absorbing state transition. Detailed 
analysis of the model reveals that the transition occur- 
ring at w c = 0.5635(1) is identical to the that observed in 
RASEP with n = 1, except that the activity now propa- 
gates symmetrically in space, looking like a familiar DP 
cluster (see Fig. [TJ. The critical exponents characterizing 
this transition also belong to the DP class; we skip the 
details of this study here. 



Summary. — Restricted exclusion processes are 
known to undergo absorbing phase transitions character- 
ized by integer exponents even though the order parame- 
ter is a fluctuating scalar with no special symmetry except 
the presence of the additional conserved density field. This 
could have been blamed for the non-DP behaviour but for 
the examples of Manna model, CTTP etc. which show 
DP behaviour even in presence of a conserved field. In 
this work wc show that when density conservation is vio- 
lated, by allowing deposition and evaporation of particles, 
the system shows an APT at a finite deposition rate w c . 
The resulting critical behaviour becomes DP, confirming 
that the conservation of density was indeed responsible for 
the non-DP behaviour in RASEP. 

The density p(w) is non-analytic at the critical point 
and we find that p — p{w c ) also undergoes a transition 
at w = w c with critical exponents identical (within error 
limits) to those of DP. 

The original transition in conserved RASEP was ob- 
served by tuning density. In the non-conserved version, 
the activity p a also show a transition w.r.t. p at p c — p{w c ) 
with exponents /?', much different from DP. We ar- 

gue that they are in fact scaled by a factor 1/Pdp- 

The one dimensional CLG, which is the symmetric ver- 
sion of a special case p = 1 of RASEP, also flows to DP 
when density conservation is broken. In a recent work, 
Fiore et. al. [20] have claimed that a modified CLG, which 
includes next nearest neighbour interaction and simulta- 
neous jump of two consecutive particles in opposite direc- 
tion but does not violate particle conservation, shows non- 



APT 


a I'll z P/v± 


DP UJ 

pa VS. W 
p VS. W 
Pa VS. p 


0.1594 0.2764 1.733 1.5802 0.2520 
0.15984(7) 0.278(1) 1.74(1) 1.58(1) 0.254(4) 
0.1596(1) 0.278(3) 1.73(1) 1.58(1) 0.256(4) 
1 6.269 



Table 1: Critical exponents of non-conserved RASEP. 

trivial critical behaviour. It still remains to see whether it 
is possible to modify the dynamics of RASEP in a differ- 
ent way, keeping density conservation, to obtain the DP 
critical behaviour. 

REFERENCES 

[1] M. Henkel, H. Hinrichsen, and S. Liibeck, Non- 
Equilibrium phase transitions, vol. 1, Springer, Berlin, 
2008; H. Hinrichsen, Adv. Phys.49, 815 (2000). 

[2] H. K. Janssen, Z. Phys. B 42, 151 (1981); P. Grassberger, 
Z. Phys. B 47, 365 (1982). 

[3] J. Marro and R. Dickman, Nonequilibnum Phase transi- 
tion in lattice models, Cambridge University Press, Cam- 
bridge, 1999. 

[4] F. D. A. A. Reis, Braz. J. Phys. 33, 501 (2003). 

[5] P. Grassberger, J. Stat. Phys. 79, 13 (1995). 

[6] M. A. Miinoz and R. Pastor-Satorras, Phys. Rev. Lett. 

90, 204101 (2003). 
[7] P. K. Mohanty and D. Dhar, Phys. Rev. Lett. 89, 104303 

(2002). 

[8] G. Odor and A. Szolnoki, Phys. Rev. E 53, 2231 (1996). 
[9] K. A. Takeuchi, M. Kuroda, H. Chate, and M. Sano, Phys. 

Rev. E 80, 051116 (2009). 
[10] I. Jensen, J. Phys. A 32, 5233 (1999). 
[11] J. W. Essam, J. Phys. A 22, 4927 (1989). 
[12] H. Takayasu and A. Yu. Tretyakov, Phys. Rev. Lett., 68, 

3060 (1992); I. Jensen, J. Phys. A 26, 3921 (1993). 
[13] M. J. de Oliveira, Phys. Rev. E 71, 016112 (2005); E. 

F. da Silva, and M. J. de Oliveira, Jnl. Phys. A : Math. 

Theor. 41, 385004 (2008). 
[14] U. Basu and P. K. Mohanty, Phys. Rev. E 79, 041143 

(2009). 

[15] M. A. Miinoz, R. Dickman, A. Vespignani and S. Zapperi, 
Phys. Rev. E 59, 6175(1999); R. Dickman, M. A. Miinoz, 
A. Vespignani and S. Zaperi, Brazilian Jnl. Phys. 30, 27 
(2000). 

[16] M. Rossi, R. Pastor-Satorras, and A. Vespignani, Phys. 

Rev. Lett. 85, 1803 (2000). 
[17] M. Basu, U. Basu, S. Bondyopadhyay, P. K. Mohanty, 

and H. Hinrichsen, Phys. Rev. Lett. (2012). 
[18] U. Basu, M. Basu, and P. K. Mohanty, submitted to Phys. 

Rev. E. 

[19] S-G. Lee and S. B. Lee, Phys. Rev. E 77, 021113 (2008). 
[20] C. E. Fiore and M. J. Oliveira, Braz. J. Phys. 36, 218 
(2006). 



